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We perform numerical simulations of the Sivers effect from single spin asymmetries in Drell-Yan 
processes on transversely polarized protons. We consider colliding antiprotons and pions at differ- 
ent kinematic conditions of interest for the future planned experiments. We conventionally name 
"framework I" the results obtained when properly accounting for the various flavor dependent po- 
larized valence contributions in the numerator of the asymmetry, and for the unpolarized nonvalence 
contribution in its denominator. We name "framework II" the results obtained when taking a suit- 
able flavor average of the valence contributions and neglecting the nonvalence ones. We compare the 
two methods, also with respect to the input parametrization of the Sivers function which is extracted 
from data with approximations sometimes intermediate between frameworks I and II. Deviations 
between the two approaches are found to be small except for dilepton masses below 3 GeV. The 
Sivers effect is used as a test case; the arguments can be generalized to other interesting azimuthal 
asymmetries in Drell-Yan processes, such as the Boer-Mulders effect. 
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I. INTRODUCTION 



High-energy collisions of (polarized) hadrons represent a testground for the theory of strong interactions, the 
quantum chromodynamics (QCD). In fact, several experiments have been performed that still await for a satisfactory 
interpretation of the resulting data. Among others, in collisions of the kind pp^ — > h^X an azimuthally 

asymmetric distribution of semi-inclusively produced hadrons h (with respect to the normal of the production plane) 
is observed when flipping the transverse spin of the target proton p or of the final hadrons h, the so-called transverse 
single-spin asymmetry (SSA). Perturbative QCD, as it can be calculated in the collincar massless approximation, 
cannot consistently accommodate these SSA sometimes very large also at high energy. 

More recently, a series of SSA measurements [fl 0, S @ m semi-inclusive lp^ — > I'nX Deep-Inelastic Scattering 
(SIDIS) has renewed the interest about the QCD spin structure of hadrons, mainly because the theoretical situation 
appears more transparent. In fact, while in hadronic collisions like pp^ — > h>^X the factorization proof is compli- 
cated by higher-twist correlators |ldj | and the power-suppressed asymmetry can be produced by several (overlapping) 
mechanisms, in SIDIS a suitable factorization theorem |lllll2| allows to clearly separate terms with different azimuthal 
dependences in the leading-twist cross section. 

The main feature of this factorization proof is the possibility of going beyond the collinear approximation, which 
opens new perspectives about the explanation of the observed SSA in terms of intrinsic transverse motion of partons 
inside hadrons, and of correlations between such intrinsic transverse momenta and transverse spin degrees of freedom. 
One of the most popular examples is the socalled Sivers effect jl3j , where an asymmetric azimuthal distribution of 
detected hadrons (with respect to the normal to the production plane) is obtained from the nonperturbative correlation 
p T xP-S T , with p T the intrinsic transverse momentum of an unpolarized parton inside a target hadron with momentum 
P and transverse polarization S T . The size of the effect is driven by a new Transverse- Momentum Dependent (TMD) 
leading- twist partonic function, the socalled Sivers function f^ T , which describes how the distribution of unpolarized 
partons is distorted by the transverse polarization of the parent hadron. Then, the extraction of f^ T would allow to 
study the orbital motion and the spatial distribution of hidden confined partons, with interesting connections with 
the problem of the proton spin sum rule and the powerful formalism of Generalized Parton Distributions |14| . 
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In single-polarized Drell-Yan processes like H x Hl -> l+l-X, there is a situation similar to SIDIS: a suitable 
factorization theorem holds 0, ^| and different asymmetric contributions can be clearly distinguished. The cross 
section must be differential in the azimuthal orientation <j) of the final lepton plane and of the hadron polarization 
cf>s with respect to the reaction plane |l6j| : at leading twist, it includes a term driven by f^~ T with the characteristic 
sin(0 — </>s) dependence. Surprisingly, there are no data for this process. At the same time, quite interestingly the 
extraction of f^ T from a Drell-Yan SSA would allow to verify its predicted sign change with respect to SIDIS [17| . a 
theorem based on general grounds which represents a formidable test of QCD universality. 

Drell-Yan measurements, with unpolarized and/or transversely polarized hadrons, are planned by several exper- 
imental collaborations (RHIC at BNL, COMPASS at CERN, PANDA and PAX at GSI, and, possibly, also future 
experiments at JPARC) . In a series of previous papers [TsL I20I l2l| , we performed numerical simulations of Drell- 
Yan SSA with transversely polarized protons using colliding protons, antiprotons, and pions, in various kinematics of 
interest for the planned experiments. In particular, we verified that the foreseen setup of RHIC and COMPASS, with a 
reasonable sample of Drell-Yan events, should allow to unambiguously extract f^ T from the correspondin g si n(qi> — </>s) 
asymmetry, as well as to clearly test its predicted sign change with respect to the SIDIS asymmetry 12(1 l21l . We 
also explored another interesting piece of the single-polarized Drell-Yan cross section [21| (see also Ref. [23] )• It is 
driven by the sin(0 + <fis) asymmetry and is related to another TMD function, the Boer-Mulders function , which 
describes the distribution of transversely polarized partons inside unpolarized hadrons. The interest in arises 
from the possibility of directly linking it to the long-standing problem of the violation of the socalled Lam- Tung 
sum rule [16( , namely the presence of an anomalous cos 2</> asymmetry in the distribution of Drell-Yan muon pairs 
in pion-induced unpolarized collisions [2^, 0, |2i| (but apparently not present in the recent data of Ref. [2(| about 
high-energy proton-deuteron collisions), which neither complicated QCD calculations at higher order, nor higher twist 
contributions, are able to justify in a consistent picture [22i 123, |2!j , and that it can alternatively be interpreted as a 
QCD vacuum effect [30j |. 

Our simulations of Drell-Yan SSA need a phenomenological input for the various TMD functions involved. Actually, 
one of our goals was to test the relation between the statistical uncertainty of simulated events and the theoretical 
uncertainty originating from different parametrizations of the same TMD function. One of the main features of these 
phenomenological analyses is the approximation of neglecting both polarized and unpolarized nonvalence partons 
(see, e.g., Ref. [3l]]), which actually amounts to effectively include their contribution in the fitting parameters of the 
valence partons for the x range considered. In the Monte Carlo code of Refs. [TsL H^. I20I l2l| . we consistently take the 
same approach, but we further conveniently make a suitable flavor average of the valence contribution which allows 
for a great simplification of formulae. We conventionally name this scheme as "framework II". Here, we consider 
also the socalled " framework I" , where we release the approximation about the flavor average and we include also 
the unpolarized nonvalence contribution. The goal is to critically discuss the two methods, also with respect to the 
approximated framework introduced by those parametrization of the Sivers function that are somewhat intermediate 
between them. More specifically, we want to identify and quantify the deviations of the results obtained within 
"framework I" from those obtained within "framework II". At a qualitative level, the origin of these deviations 
can be easily identified a priori, while from the quantitative point of view there are quite distinct situations that 
need to be separately analyzed. We will systematically adopt the Sivers effect as our test case, because the relative 
abundance of these SSA data allows to already build realistic parametrizations of f^~ T , some of which are discussed in 
Sec. [HJ However, the argument can be generalized, in principle, also to the azimuthal asymmetry generated by the 
Boer-Mulders effect or by the violation of the Lam- Tung sum rule. 

The paper is organized as follows. In Sec.^ the general formalism, the considered phenomenological parametriza- 
tions of the Sivers function, and the main approximations leading to the definition of "framework I" and "framework 
II" , are discussed. In Sec. 11111 simulations for SSA within the two frameworks are compared for different Drell-Yan 
collisions on transversely polarized proton targets in several kinematic conditions of interest. Finally, in Sec. II VI some 
conclusions are drawn. 




FIG. 1: The Collins-Soper frame. 
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II. GENERAL FORMALISM AND APPROXIMATIONS 



In a Drell-Yan process, an antilepton and a lepton with individual momenta k\ and ki are produced from the 
collision of two hadrons with momentum P iy mass Mj, and spin Si, with i = 1,2. The center-of-mass (cm.) square 
energy available is s = (Pi + P2) 2 and the invariant mass of the final pair is given by the time-like momentum transfer 
q 2 = M 2 — (ki + k2) 2 ■ If M 2 , s — > 00, while keeping the ratio < r = M 2 /s < 1 limited, the factorized elementary 
mechanism proceeds through the annihilation of a parton and an antiparton with momenta p\ and P2 , respectively, 
into a virtual photon with time-like momentum q 2 . If P^ and P 2 _ are the dominant light-cone components of hadron 
momenta in this regime, then the partons are approximately collinear with the parent hadrons and carry the light-cone 
momentum fractions < x\ — p\ / P^ , X2 = P 2 /P 2 ~ < 1) with q + = p^ , q~ = p^ by momentum conservation |16| . 
The transverse components Pi T oipi with respect to the direction defined by Pi(i = 1,2), are constrained again by the 
momentum conservation q T — p\ T +P2t, where q T is the transverse momentum of the final lepton pair. If q T ^ the 
annihilation direction is not known. Hence, it is convenient to select the socalled Collins-Soper frame [32j described 
in Fig.^ The final lepton pair is detected in the solid angle (0,<fi), where, in particular, <f> (and all other azimuthal 
angles) is measured in a plane perpendicular to the indicated lepton plane but containing h = q T /\q T \. 

By neglecting terms ~ 1/M 2 , with M the largest mass of the initial hadrons, the expression of the leading- twist 
differential cross section for the H1P1I — > X process can be written as [lfij 
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where a is the fine structure constant, dVL — sin 9d8d(j), e q is the charge of the parton with flavor q, </>g 2 is the azimuthal 
angle of the transverse polarization vector of the hadron H\ , and 
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2~y + lr) = 4(1 + cos 2 1 



B{y) =y(l-y) = - sin 2 



(2) 



The TMD functions /' (H), hf q (H), describe the distributions of unpolarized and transversely polarized partons in 
an unpolarized hadron H, respectively, while f^~ T q (H^) and h\(H^) have a similar interpretation but for transversely 
polarized hadrons H* . The convolutions are defined as 



TMDf (ifi) TMD|(^ 



dpi T dp 2T S {pit + P2t - q T ) 

X [TMT) 1 (x 1 ,p lT -q/H 1 )TMr>2{x2,P2T;q/H^ ) ) + (q^q) . (>',) 
The Monte Carlo events have been generated by the following cross section fl8j : 
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(4) 



where the event distribution is driven by the elementary unpolarized annihilation, whose transition amplitude T has 
been highlighted. In Eq. Q, we assume a factorized transverse- momentum dependence in each TMD such as to break 
the convolution T , leading to 
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where q T = \q T \. The function A is parametrized and normalized as in Ref. |25| . where high-energy Drell-Yan tt — p 
collisions were considered. The average transverse momentum turns out to be (q T ) > 1 GeV/c (see also the more 
recent Ref. 01), which effectively reproduces the influence of sizable QCD corrections beyond the parton model 
picture of Eq. |Q. It is well known [34J that such corrections induce also large K factors and an M scale dependence 
in parton distributions, determining their evolution. As in our previous works |l8l Il9l l2dL l2l|. we conventionally 
assume in Eq. J3J that K = 2.5, but we stress that in a single- spin asymmetry the corrections to the cross sections 
in the numerator and in the denominator should compensate each other, as it turns out to actually happen at RHIC 
cm. square energies [35|. Since the range of M values here explored is close to the one of Ref. [25|, where the 
parametrization of A,F, and c,; in Eq. ffl, was deduced assuming Af-independent parton distributions, we keep our 
same previous approach [H, l23> Sll an d use 
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(6) 



where the unpolarized distribution ff(x) for various flavors q — u,d, s, is taken again from Ref. |2q . 

The whole solid angle (0, <fi) of the final lepton pair in the Collins-Soper frame is randomly distributed in each 
variable. The explicit form for sorting it in the Monte-Carlo is [la, 120, l2Jj 
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If quarks were massless, the virtual photon would be only transversely polarized and the angular dependence would 
be described by the functions c\ = 1 and Si = 1 + cos 2 9. Violations of such azimuthal symmetry induced by the 
function c 2 = ^ are due to the longitudinal polarization of the virtual photon and to the fact that quarks have an 
intrinsic transverse momentum distribution, leading to the explicit dependence of v upon q T and to the violation of 
the socalled Lam- Tung sum rule [2^ . I2~H l25l | . QCD corrections influence v, which in principle depends also on M 2 [2^|. 
Azimuthal cos 2(f) asymmetries were simulated in Ref. [T^ | using the simple parametrization of Ref. [Tfil | and testing it 
against the previous measurement of Ref. [2^, 13 125| . 
The next term in Eq. J7J) describes the Sivers effect [T^ |: 



E 9 



C.3 = Csiv(q T ,Xl,X 2 ) 



h ■ p 2j 



M 2 



E 9 e%F[fi(.Zi,Pi T ) ff{x 2 ,p 2T ) 



(8) 



while the last one contains the asymmetry induced by the socalled Boer-Mulders effect [l6l ]: 
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In the following, we will consider the asymmetry generated only by csiv in Eq. (jSJ), because there is a sufficient 
amount of available data to support the construction of realistic parametrizations for the Sivers function f^. However, 
our arguments can be easily generalized also to the Boer-Mulders cbm term. We will come back on the v coefficient 
at the end of next Section. 

For sake of consistency, the denominator of Eq. JSJ) is approximated by the same |T| 2 of Eq. (JSJ. As for the 
numerator, we first simulate the Sivers effect using the parametrization of Ref. |3l| , 
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where M 2 is the mass of the polarized proton, p T = \p T \, and (p 2 ) = 0.25 (GeV/c) 2 is deduced by assuming a Gaussian 
ansatz for the p T dependence of fx in order to reproduce the azimuthal angular dependence of the SIDIS unpolarizcd 
cross section (Cahn effect). Flavor-dependent normalization and parameters in the x dependence are fitted to SIDIS 
SSA data using only the two flavors q — u,d and neglecting the (small) contribution of antiquarks (see Refs. [2TL l3l| ^ . 
Following Ref. 0] and including a sign change of when plugging it in the Drell-Yan cross section 0], we get 
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where, for brevity, Ug iv (x) represents the contribution of flavor u to the x dependence of the Sivers function 
parametrized as in Eq. (|1C)|> (and similarly for flavor d). 

As an alternative choice, we adopt the new parametrization described in Ref. |20l l2l| . It is inspired to the one of 
Ref. [3(j , whose x dependence is retained but a different flavor-dependent normalization and an explicit p T dependence 
are introduced. The latter is bound to the shape of the recent RHIC data on pp^ — > ttX at i/s — 200 GeV |37j . 
where large persisting asymmetries are found that could be partly due to the leading-twist Sivers mechanism. The 
expression adopted is 
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where po = 2 GeV/c, and N u = —Nd = 0.7. The sign, positive for u quarks and negative for the d ones, already takes 
into account the predicted sign change of fx T from Drell-Yan to SIDIS ^3 • 
Along the same previous lines, we get 
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where now Ug iv (x) indicates the contribution of flavor u to the x dependence of the Sivers function parametrized as 
in Eq. (Pf. 

In the following, we will refer to "framework I" as to the (1 + cos 2 9) sin(0 — 4>s 2 ) angular asymmetry generated in 
Eq. CJ by the coefficients JTTJ or (0 . 

The coefficients Cg iv and Cg iv can be further approximated with a procedure that here we will conventionally indicate 



as "framework II". Again, following the lines described in Refs. |20j,|21j we obtain 
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The coefficients n q (q = u, d) include the contribution of the quark charge, of the normalization of the parton 
distributions, as well as of the statistical weight of the considered flavor. In fact, the approximation is based on the 
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idea that each term in both flavor sums in the numerator and denominator can be replaced by a "flavor- averaged" 
one, and the resulting simplified ratio is then properly weighted. For example, in a pp^ collision where only valence 
contributions are considered and the u(x) distribution in p is normalized to 2 (and similarly for the u one in p), we 
easily get a statistical ratio 16:1 of the wuJ annihilations over the dd) ones, such that n u — 16/17, rid = 1/17. 

The underlying idea in "framework II" is that there is no strong flavor dependence in the sums appearing in 
Eqs. (|ll|l and (|13|) . and that it is possible to neglect the contribution from sea (anti)quarks. The latter feature is 
included by construction in the parametrizations of the Sivers functions, in agreement with the common belief about 
the behaviour of SSA at low x. It is far less obvious that the approximation can be safely carried on also in the 

A/ B 

denominator of the asymmetries c Siv , where the unpolarized distributions f\ are involved. In the next Section, wc 
will numerically simulate n^p^ and pp^ Drell-Yan collisions to explore the effect of neglecting such terms. 



III. MONTE CARLO SIMULATION AND DISCUSSION OF RESULTS 



In this Section, we present results for the Monte Carlo simulation of the Sivers effect in several Drell-Yan events 
using transversely polarized proton targets by adopting either " framework I" or " framework II" , as they are described 
in the previous Section. The goal is to explore the sensitivity of the results when the contribution of unpolarized sea 
(anti) quarks is neglected in the denominator of the SSA, as it is usually done when parametrizing the Sivers function 
from experimental data. 

In the Monte Carlo, events are simulated by the cross section with "framework II" , namely using Eqs. I|14l) or l|15|) 
according to the input parametrization selected for the Sivers function |20ll2lLl3l| ; the unpolarized parton distributions 
are parametrized following Ref. 2jj, as explained in the previous Section when discussing Eqs. (JSJ and ijfjj. Events 
for the Sivers effect are then rejected/accepted by constructing the spin asymmetry A! = (U' — D')/(U' + D'), where 
U' {D') represents events with positive (negative) values of sin(</> — 4>Si) m Eq. (0. Similarly, the spin asymmetry 
A = (U — D)/(U + D) is obtained by rejecting/accepting events using the "framework I", namely using Eqs. 111|) 
or (|13l) . for positive (U) or negative (D) values of sin(</> — <fis 2 )- Data are accumulated only in the x 2 bins of the 
polarized proton, i.e. they are summed over in the x\ bins for the hadronic beam, in the transverse momentum q T 
of the lepton pair and in their zenithal orientation 0. In the following, plots will compare A and A' for different 
Drell-Yan processes and kinematical conditions, showing only the positive values of the asymmetries (the negative 
ones look equally distributed in a symmetric way). 

All events refer to collisions at s = 100 GeV 2 , which can be explored either at GSI in the socalled collider mode, 
or at COMPASS with fixed targets. The lepton pair invariant mass is constrained in the ranges 4 < M < 6 GeV 
and 1.5 < M < 2.5 GeV in order to avoid overlaps with the resonances of the cc and bb quarkonium systems, while 
exploring at the same time different (x 2 ) regions which stress the role of sea (anti)quarks. The transverse momentum 
of the lepton pair is constrained in the range 1 < q T < 3 GeV/c in order to avoid a strong dilution of the SSA because 
of the rapid decrease of the distributions (|10fl and (|12l) at larger q T . Moreover, the resulting (q T ) ~ 1.8 GeV/c is in 
fair agreement with the one experimentally explored at RHIC [3?| . 

The considered Drell-Yan collisions involve transversely polarized protons and different hadronic probes: antiprotons 
(p) and pions (ir~ and tt + ). The statistical sample is made of 7 000 events except for the ir + probe, where we have 
used 20 000 events because the Monte Carlo indicates that the cross section involving ir + is statistically disfavoured 
by approximately a factor 1 /4 [3£| . 

As a sort of "warm up", we first consider oversimplified cases where the simulated comparison between A and 
A' can be confronted with a priori known results. In Fig. |2 we show the scatter plot of A versus A' for the 
pp 1 ' — > fi + /j,~X process with 1.5 < M < 2.5 GeV where only (polarized) valence contributions are considered and 
usiv(x) = u(x), dsiv(x) = d(x). From Eq. I|14|) . it is evident that both "framework I" and "framework II" must give 
the same result, namely 

el u{x 1 )us lv {x 2 ) + e 2 d d(xi)d S iv{x 2 ) _ 



e 2 u u(x\) u(x 2 ) + e 2 d d(x\) d(x 2 ) 
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n d = l, (16) 



because any flavor dependence has been switched off and for "framework I" the asymmetry A is not diluted by the 
contribution of unpolarized sea (anti)quarks showing up in the denominator of the first equality in Eq. l|14fl itself (the 
socalled "sea dilution effect"). From the figure we get the consistent picture that the SSA calculated with the two 
methods are statistically equal. We have checked that the same result holds also for 7r ± probes, as it is obvious. 
In Fig. |3 the scatter plot is shown for the ir~p^ — > /i + pT X process with 4 < M < 6 GeV where still usi V (x) = 
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FIG. 2: Scatter plot for single-spin asymmetries A and A', calculated with "framework I" and "framework II" (see text), 
produced by 7000 events of the Sivers effect in the pp^ — > /i + n~ X process with muon invariant mass 1.5 < M < 2.5 GeV 
and transverse momentum 1 < q T < 3 GeV/c, where only (polarized) valence contributions are considered and usiv(x) = 
u{x), dsiv(x) = d(x). 
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FIG. 3: Scatter plot in the same conditions as previous figure, but for the tv p' —* fi + fi X process with muon invariant mass 
4 < M < 6 GeV and including the contribution of unpolarized sea (anti)quarks. 




u(x), dsivix) = d(x), but the contribution of the unpolarized sea (anti)quarks is included. Therefore, we have 
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[el u(xi)u(x 2 ) + e 2 d d{x\) d(x 2 ) + e 2 s(xi) s(x 2 )] + (1 *-> 2) 



u(xi)u(x 2 ) + e d d(x 1 )d(x 2 ) + e 2 s s(xi) s(x 2 ) + e 2 s{x\) s(x 2 ) 



e\ u{x\) u(x 2 ) + e d d{x{) d{x 2 ) 



< 1 



c ¥iv (X 



USiv(x 2 ) dsiv(x 2 ) 
n u ; — ; h n d 



u(x 2 ) 



d(x 2 ) 



rid = 1 



(17) 
(18) 



Evidently, the difference between the two approaches stems from the unpolarized sea-(anti)quark contribution, both 
in the numerator and in the denominator of c^, and it is responsible for the (limited) spreading of some events in 
the scatter plot. 

In Fig. 01 the same situation is shown for muon invariant masses 1.5 < M < 2.5 GeV. At the given s = 100 GeV 2 , 
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FIG. 4: Scatter plot in the same conditions as previous figure, but for muon invariant mass 1.5 < M < 2.5 GeV. 



lower M implies testing lower Xi values, typically below 0.05-0.1, where sea effects start dominating over the valence 
contribution. Consequently, the "sea dilution effect" becomes stronger and the marked difference between "framework 
I" and "framework II" produces the scattering of events observed in the plot. In summary, the two approaches become 
equivalent when the polarized distributions weakly depend on flavor (or, equivalently, one flavor dominates the flavor 
sum in the numerator of the SSA), and the contribution of sea (anti)quarks becomes negligible, i.e. for not too low 
(x) and, consequently, M (assuming that also transversely polarized sea (anti)quarks can be neglected, as well). 

a 
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0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 
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FIG. 5: Scatter plot in the same conditions as previous figure, but for the u p' — > fJ. fi X process. 




In Fig.[5l the same kinematical situation of the previous figure is considered for tt + probes. Equations ijTY|) and itl^l) 
formally still hold. However, while the latter remains obviously unchanged, the former can be usefully compared in the 
two cases by adopting some simple approximations. In fact, if we neglect the product of two sea-quark distributions and 



we assume isospin symmetry among the valence distributions of the n probes (i.e., (x) 
and d n + (x) — d v - (x) = d^ix) = d(ir)), we get for the tt~ probe 



- (x) = u Tr (x) = u(tt) 



-*Siv 7T 



d(TT) d(p) 



u(w) u(p) + d(ir ) d(p) 



(19) 



and for the 7r + probe 



C SlVTT+ <X 



4u(ir) u{p) 



d(ir)d(p) +8u(n+)u(p) 



(20) 
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For sake of simplicity, we further assume the flavor independence of the sea contributions, namely d(ir~) — u(tt + ) = 
Qseai^) "C u(tt), d(w) and u(p) = d(p) = q sea {p) <C u(j>), d(p). It is also reasonable to consider u(ir) > d(ir) and 
u(p) > d(p). From these inequalities it follows that 

d(n) q S ea(p) < 4m (tt) g sea (p) 

u(tt)u(p) + q sea (ir) d{p) d(n) d(p) + 8q sea {n) u(p) 

namely c^ ilJ7r _ is closer to c ! g iv — 1 than Cg fojr+ . This fact is responsible for the plot in Fig. [5] showing much more 
scattered events than in Fig. 01 Given the lowest M range here explored, the situation displayed in Fig. represents 
a sort of ideal upper limit to the potential discrepancy between the two methods for calculating the Drell-Yan SSA. 

a 
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0.02 

°0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 

A' 



FIG. 6: Scatter plot in the same conditions as Fig. 0] but using the parametrization of the Sivers function from Ref. |3lj . 

We now turn to the comparison between " framework I" and " framework II" with realistic parametrizations of the 
Sivers function and including the " sea dilution effect" . In Fig. HO the scatter plot of A versus A' is shown for the 
7r~p^ — » /!+ pr X process with 1.5 < M < 2.5 GeV using the Sivers function of Ref. [U, i.e. using c§ iv from Eqs. 
and lj*L4*)l . The discrepancy is not large, if we consider the potentially dangerous low M range. 
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FIG. 7: Scatter plot in the same conditions as previous figure, but using the parametrization of the Sivers function from 

Ref. mm and for muon invariant mass 4 < M < 6 GeV. 

In Fig. the same situation of the previous case is simulated at 4 < M < 6 GeV using the Sivers function from 
Ref. |20ll21| . Here, the agreement between the two methods is much more evident also because the higher M range 
considered corresponds to a safer valence domain in x%, approximately around 0.2-0.4. 

Next, in Fig. [S] we reconsider the situation of the previous figure but for the pp^ — > /i + /i~ X process. At the level 
of leading valence contributions to the asymmetry, the 7r~p^ and the pp^ collisions are equivalent because they are 




1 j>0 

-S, " i 1 1 i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 i 1 1 1 1 1 1 1 1 1 1 1 1 1 1, 
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FIG. 8: Scatter plot in the same conditions as previous figure, but for the pp T — * /i + fi X process. 

both dominated by the uu^ annihilation. Hence, the persisting close similarity of A and A 1 indicates that the main 
origin of the discrepancy between " framework I" and " framework II" comes from the " sea dilution effect" , namely 
from including or neglecting the contribution of unpolarized sea (anti)quarks in the denominator of the SSA. 

a 

0.5 r 1 

0.45 ~- 

0.4 E- 
0.35 t 

0.3 E- 
0.25 

0.2 
0.15 
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0.05 | .tflgtiPWPf 

n— w^ - **^^ i , i 

0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 
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FIG. 9: Scatter plot in the same conditions as Fig. [5] but for the 7r + p T — > fi + fj, X process. 

Last, in Fig.[§]we show the scatter plot of A versus A' in the same conditions as in Fig.HJbut for the tt + p^ — ► pi + fi~X 
process. The large discrepancies confirm the findings about Eq. I|21|) that were displayed in Fig. [S] using a simpler 
naive parametrization of the Sivers function: with the ir + probe the "sea dilution effect" is emphasized and the 
approximations introduced in "framework II" with Eq. (|14|) are somewhat not justified. The effect is also emphasized 
by the low M range considered. 

In summary, when the spin asymmetries A and A' are simulated using realistic parametrizations of the Sivers 
function, which are all based on the dominance of the transversely polarized u' valence distribution in the transversely 
polarized p^ parent hadron, their results are very close provided that the effect of the unpolarized sea (anti)quarks 
can be neglected. This condition is fulfilled by the n~ and p probes at not too low muon invariant masses M, which 
correspond to the safe valence domain in x. In these cases, the originating schemes "framework I" and "framework 
II" can be considered equivalent. 

However, the last statement can be misleading, because the realistic parametrizations of the Sivers function are 
obtained in the "framework II", namely by building the SSA by neglecting the unpolarized sea (anti)quarks. It 
means that any contribution from the sea, which is intrinsically contained in the experimental measurement of the 
asymmetry, is effectively reproduced in the parametrization, particularly in the behaviour at low x. For example, this 
" hidden sea effect" is contained in the a q parameter of the x aq factor in Eq. I|ll|) , which is determined from a fit of 
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experimental SSA based on valence q = u,d degrees of freedom only |3l|. Hence, it can be over- or under-estimated 
when the related Sivers function is plugged into expressions of SSA constructed with "framework I" or "framework 
II". 

Actually, the emerging issue here is that "framework I", which looks more correct because it properly includes the 
unpolarized sea contribution in the denominator of the SSA, could produce a sort of double counting when employing 
a Sivers function in the numerator that is parametrized with valence quarks only. Consequently, it is not obvious that 
the approximations adopted in "framework II" (and systematically used in the analysis of several Drell-Yan SSA in 
Refs. |!8lll9ll2Ctl2l| ) are crude ones; rather, they look like the most appropriate framework for using the parametrized 
Sivers functions available in the literature. However, this last statement should be generalized with some care. In 
fact, the large discrepancies observed in the results obtained with "framework I" and "framework II" using ir + probes, 
suggest that each physics case should be separately considered. 

Finally, as anticipated in Sec. [Til we reconsider the v term in Eq. l(7)l. which is most likely related to the violation 
of the Lam- Tung sum rule |23l l24l l25j| . The corresponding azimuthal cos 2cf> asymmetry in the unpolarized Drell- 
Yan cross section was simulated in Ref Hall using the simple parametrization of Ref. and testing it against the 
previous measurements of Refs. [23l E3. l25l ] . The Boer-Mulders function h^(x) was parametrized in Ref. 0] exactly 
in the "framework II" but with no flavor dependence, because of few available data at that time. The latter span 
only the region x > 0.3, while for x < 0.3 the v function was assumed to be independent from x. Hence, the present 
parametrization of h\ does not effectivel y in clude any "sea dilution effect" . The situation could be improved by either 
refitting the Drell-Yan data of Refs. |23l 124. |25| using "framework I" and an ^-independent v function for x < 0.3, or 
using "framework II" and introducing a specific x a power law to extrapolate the behaviour at low x. Moreover, new 
data have been recently released for Drell-Yan muon pairs produced in high-energy proton-deuteron collisions |26| . 
that show no evidence of a cos 20 asymmetry at very low x. In any case, the few available data do not allow yet a 
full flavor-dependent analysis and any conclusion is, consequently, premature. 



IV. CONCLUSIONS 

In this paper, we performed numerical simulations of the socalled Sivers effect |l3j . as it can be isolated in single 
spin asymmetries (SSA) of the distribution of Drell-Yan muon pairs produced from collisions of transversely polarized 
protons and different hadronic probes. Several measurements of such SSA are planned by experimental collaborations 
(RHIC at BNL, COMPASS at CERN, PANDA and PAX at GSI, and, possibly, also future experiments at JPARC). 
The goal is the extraction of the socalled Sivers function /j^,, a leading- twist parton distribution that can give direct 
insight into the orbital motion and the spatial distribution of hidden confined partons, with interesting connections 
with the problem of the proton spin sum rule and the powerful formalism of Generalized Parton Distributions . 
In particular, it would be extremely important to verify the peculiar universality properties of such parton density, 
namely its predicted sign change with respect to the same f^ T as it is parametrized to reproduce the SSA in SIDIS 
measuremen ts [ 20U2lj| . This theorem is based on very general assumptions and it represents a formidable test of QCD 
universality |17|. 

One of the main features of the SIDIS phenomenological parametrizations is the approximation of neglecting both 
polarized and unpolarized nonvalence partons (see, e.g., Ref. H3|), which actually amounts to effectively include their 
contribution in the fitting parameters of the valence partons for the x range considered. We conventionally name this 
the "hidden sea effect". In a series of previous papers 0, EJUS El > we performed numerical simulations of Drell- 
Yan SSA with transversely polarized protons using colliding protons, antiprotons, and pions, in various kinematics of 
interest for the planned experiments. In our Monte Carlo code, the Sivers effect was consistently simulated within 
the same approach, but we further conveniently made a suitable flavor average of the valence contribution which 
allows for a great simplification of formulae. We conventionally name this scheme as "framework II". Here, we 
consider also the socalled " framework I" , where we release the approximation about the flavor average and we include 
also the unpolarized nonvalence contribution, which shows up mainly in the denominator of the SSA; as such, we 
conventionally refer to the " sea dilution effect" . 

We have explored the deviations of "framework II" from the more appropriate "framework I" in the simulation of 
the Sivers effect for the Hp^ — > fi + fi~X process at s = 100 GeV 2 with H = p,ir~ ,tt + and different muon invariant 
masses. It turns out that the two approaches become approximately equivalent when the polarized distributions 
weakly depend on flavor (or, equivalently, one flavor dominates the flavor sum in the numerator of the SSA), and 
for sufficiently high muon invariant masses (typically, 4 < M < 6 GeV), which correspond to (x) in the safe valence 
domain where the contribution of sea (anti)quarks becomes negligible (assuming that also transversely polarized sea 
(anti)quarks can be neglected, as well). These conditions can be fulfilled when using the n~ and p probes; on the 
contrary, the valence structure of ir + emphasizes the role of the "sea dilution effect" and leads to larger discrepancies. 

However, we must remark that the parametrizations of the Sivers function are basically obtained using "framework 
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II", since all sea (anti)quarks are neglected. Hence, the "hidden sea effect" contained in the values of the parameters 
can be over- or under-estimated when the related Sivers function is plugged into expressions of SSA constructed with 
"framework I" or "framework II". Actually, the "framework I" appears even less appropriate, since it could lead to a 
double counting of the contribution of unpolarized sea (anti)quarks; in our jargon, the "sea dilution effect" induced 
by the denominator of the SSA would describe the same mechanisms as the "hidden sea effect" contained in the 
fitting parameters of the valence quark distributions. Correspondingly, the approximations adopted in "framework 
II" (and systematically used in the analysis of Refs. |18L figl |20L l21| ^ look like the most appropriate approach for using 
the parametrized Sivers functions presently available in the literature. However, this statement should be taken with 
some care, because there are cases like the ir + p^ collisions where the "framework II" largely deviates from "framework 
I" in any kinematics and seems not well justified. 

In our work, the Sivers effect has been used as a test case since the abundance of data allows for a flavor-dependent 
analysis. In principle, the arguments can be generalized to other interesting azimuthal asymmetries in Drcll-Yan 
processes, such as the Boer-Mulders effect or the violation of the Lam- Tung sum rule. But in the former case, there 
are no experimental data, while in the latter the few ones available [13)13 E3- do not permit to discriminate the 
contributions of each flavor and prevent from coming to definite conclusions. 
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